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■ ■ - I. INTRODUCTION 

00 ■ 
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' Using a very simple model with asymptotic freedom (or limit cycle) and bound states, this lecture explains a 
fSJ . renormalization group (RG) procedure for Hamiltonians, including the result that an effective interaction strength 
' grows at low energies. Instead of integrating out high-energy degrees of freedom, the procedure explained here is carried 
, out using similarity transformations. The magnitude of increase of the interaction strength (a coupling constant in 
' an effective Hamiltonian) at low energies is related in the model to the existence of a bound state and a degree to 
, which the interaction responsible for the effect of binding is included in the generator of the similarity transformation. 
^Sj ■ Most briefly, the more interaction in the generator the less increase in the coupling constant. In addition, since a tiny 
' and quite generic modification of the simple model changes the asymptotically free behavior into a limit cycle of a 
' _^ '. large period, the model shows that an apparently asymptotically free behavior in a considered range of scales may 
' actually be a part of a cycle with a period much larger than the considered range. In the limit cycle case, the model 
^ , illustrates how a large set of bound states can influence behavior of the effective coupling constant at low energies, 
^ ' also depending on the generator used for the similarity transformation. Key references for this lecture are [l[ and 
] How the renormalization group findings concerning bound states may be related to the AdS/CFT correspondence [3|, 
0], including the RG interpretation proposed by Polyakov [1, @], is not yet established. On the other hand, the RG 
I I procedure discussed here is particularly useful in the front form of Hamiltonian dynamics [3, ■ This form of dynamics 
^ I is precisely the one in which recent holo grap hic models have been proposed for calculating masses and wave functions 
(yQ ' of bound states of quarks and gluons [3, [l^l ■ 

I/"} : Section[TT]defines the model Hamiltonian. Subsequent sections introduce relevant concepts on the basis of elementary 
C^l ' calculations in the model. Ultraviolet divergences are identified in Section IIIII Section IIVI reviews a standard RG 
, procedure. The review includes: derivation of the result that a coupling constant increases at low energies (without 

■ explicit discussion of scattering processes or Green's functions) ; demonstration of asymptotically free behavior of the 

■ coupling constant; explanation of how counterterms are calculated; and identification of difhculties associated with the 
00 i increase of the coupling constant at low energies. Section |V] introduces the similarity RG (SRG) procedure. A simple 

' class of generators of similarity transformations is introduced in Section IV Al Solutions for effective Hamiltonians in 
^ , the SRG procedure with different generators, are described in Section IVTl The case of asymptotic freedom is discussed 
in Section IVI Al and the case of limit cycle, related to asymptotic freedom, is discussed in Section IVI Bl Section lylll 
^ , concludes the lecture by a summary of the connection found in the model between the increase of effective interaction 
' strength at small energies and existence of bound states. 

A set of brief appendices is added in order to point out analogies between the simple model and quantum field 
theory (QFT), such as QCD. The appendices discuss: singular, (5-function potentials in effective Schrodinger equations 
(Section I VIII A[) . theory of effective particles that aims at explaining the origin of the constituent quark model and 
binding of partons (Section IVIIIB[) . and some questions concerning a possible connection between the SRG procedure 
and AdS/CFT correspondence (Section IVlirC)) . 
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II. MODEL 



In order to understand the mechanism that causes an asymptotically free interaction to rapidly increase at small 
energies, a theoretical feature often suggested to be associated with confinement, we need a precise definition of what 
is meant by the effective coupling constant. In theories as complex as QCD, the concept of effective coupling constant 
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at low energies is not simple. Therefore, it seems appropriate to first define the concept using a simple model and see 
what happens there before one plunges into the complexity of QCD. 

The well-known asymptotically free coupling constant depends on the energy scale logarithmically (it is proportional 
to the inverse of the logarithm of the energy). The energy is actually a kinetic one, not involving interactions in any 
significant way, and the scale should be conceptually associated with momentum, a kincmatical variable of quantum 
mechanics rather than a dynamical one. Thus, when constructing a simple model Hamiltonian in the generic form of 

H = Ho + Hi, (1) 

one considers Hq a free (kinetic) energy and Hj an interaction. An alternative way of thinking is that all eigenvalues 
and eigenstates of Hq are known and provide a basis for describing what happens due to Hj. In both ways, for studies 
of logarithmic effects of asymptotic freedom, it is useful to assume that the spectrum of Hq is not degenerate and has 
the form 

HqIu) = En\n), En^tih", 6> 1, (2) 

where /i is a unit of energy. For convenience, is set to 1 and omitted from further consideration (all numbers that 
refer to energy need to be multiplied by ^ in order to obtain the actual quantity). The utility of using powers of 
h is that successive momentum scales, or eigenvalues of Hq, are separated by a constant on a logarithmic scale, the 
constant being Info. Therefore, a small number of states, just one per momentum (or energy measured by Hq) scale, 
will be sufficient to track logarithmic effects associated with asymptotic freedom and, at the same time, control effects 
of binding. Once the eigenstates of Hq are normalized, (m|n) = 5mm the matrix elements of Hq in the model are 

{m\HQ\n) = EmSmn- (3) 

The interaction Hamiltonian Hi is defined in the model by its matrix elements, using assumption that these elements 
should be factorized, i.e., Hjmn — HmHn, for the purpose of having very simple exact solutions to the eigenvalue 
problem for H. Thus, the factor i7„ should have dimension of square root of energy. The simplest possibility is 
Hn ~ \fEn- The proportionality is reduced to a dimensionless number, and one can write 

(MUM) = -a^^E^^^fE^^, (4) 

where g determines the strength of the interaction. It is called coupling constant, in analogy with the standard 
nomenclature in QFT. The negative sign is chosen so that there exists a bound state for a sufficiently large positive 
g. By definition, the bound state corresponds to a negative eigenvalue of H . 
In summary, matrix elements of the model Hamiltonian are 

Hran = {m\H\n) = E,nSmn - 9 \/eZ\/E^ ■ (5) 

Note that it will have to be clarified what the words "sufficiently large g" are supposed to mean, because the model 
Hamiltonian produces divergences (infinities) no matter how small the number g is. 



III. ULTRAVIOLET DIVERGENCES 

Consider first very small numbers g and the intuition that eigenvalues of H should be nearly equal to the eigenvalues 
of Hq. For example, first-order perturbation theory produces a correction AEm^ to the energy E„i of the form 

AS« = {m\Hi\m) ^-gEm- (6) 

When g is small, the correction is small; it is just a fraction g of the energy being corrected. (This situation resembles 
what happens in the Schrodinger equation for atoms or positronium when one derives corrections to the Coulomb 
potential from QED and calculates their influence on energy levels in first-order perturbation theory, see Appendix, 

Section EmAl) 

Consider now the second-order correction to the same energy 

AEi^) = y MEiMl . g^Em y . (7) 

Em — Ek f—^ Em — Ek 

Since Ek = b'^, a term number k contributes 1/(6™^'^ — 1) in the sum. Thus, terms with k < m contribute a finite 
sum of a nearly geometric series (with quotient 1/6), which is not sensitive to any lower bound on k in the sum, say 
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M, if ^ 6™. However, a lower bound is needed to define a Hamiltonian whose matrix has a finite size. In order 
to focus attention on Em ~ 1, one can set the lower bound on A; to be a large negative integer M. As a result, all 
eigenvalues of Hq included in the model are from now on assumed to satisfy the condition Ek > b^^ , and 6*^ ^ 1. 

In contrast, terms with k ^ m contribute each —1. For 6 S> 1, the second-order correction is proportional to the 
number of basis states with energies (eigenvalues of Hq) greater than Em- In order to obtain a finite answer, one must 
limit this number. For example, one can limit from above the range of summation over k by certain large positive 
integer N. Then, one has to understand what happens when N is very large. Imposing a limit such as N is called 
ultraviolet regularization. 

Including the ultraviolet regularization, the result of Eq. ([7]) is approximated by the formula (the larger b, the 
better the approximation) 

AeI^'> - -g^E,n{N-m). (8) 

This formula shows that the correction tends to — oo when one attempts to send the ultraviolet cutoff on energies, 
A = b^ , to infinity. Since N = In A/ In 6, the second-order correction is also approximately given by the formula 

This result explains why the correction is called logarithmically divergent in the ultraviolet. A direct comparison 
of Eqs. ([HI) and ^ illustrates that the ultraviolet logarithmic divergence results from all different energy scales 
contributing equally to the eigenvalues. This conclusion is not limited to perturbation theory. 
Consider the eigenvalue problem for the matrix [Hmn] , 

TV 



Hrnnlpn = Elpm , (10) 




n=M 

whose subscripts are limited after the regularization, M < m,n < N . Because the interaction is factorized, one has 

V'm = " V^V-n (11) 

and the sum in Eq. (Ilip is just a number (does not depend on m), say c. Substituting this solution for the wave 
fimction into Eq. (fTO| . one obtains 



(12) 

E„i — E 



As far as the high-energy part is concerned, the sum in this non-perturbative eigenvalue condition is of the same type 
as in the perturbative Eq. ([7]), the only change being that Em with one particular value of m is now replaced by the 
unknown eigenvalue E. This does not change the fact that the sum can be approximated by — In (A/ji?!). The result 
is that the eigenvalue condition has no meaning for finite E when A — > oo. 

More specifically, consider a possibility that a bound state exists, with a negative eigenvalue E — —Eb- Replace 
the sum over n by an integral with measure dn = dE„ / {En In b), just to get an idea what happens. One has 

q dE 

or 



In b Jf,M Eb+E 



Eb - . (15) 



e 9 — 1 

When one takes the limit of A — s- oo for fixed g, the binding energy diverges linearly with A. 
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One can also observe that the square of the matrix Hi with large N and M is equal to Hj , which means that 
it diverges in the ultraviolet limit of A — > oo for fixed g. This means that all powers of the entire H are ultraviolet 
divergent. In particular, the evolution operator [/(t, 0) = e^'^* does not exist in this limit. 

Note that the divergences in Eqs. ([5]), and p^ . result from the diverging number of degrees of freedom. Similar 
problems may occur in classical statistical systems with a huge number of degrees of freedom if all these degrees of 
freedom contribute significantly to observables [lll |. 

IV. STANDARD RENORMALIZATION GROUP PROCEDURE 

Given the divergences, the task now is to figure out whether to keep or discard the model that produces such 
diverging results when A'' grows. Of course, one can keep N fixed and try to describe physics (in the simple model, 
"physics" amounts to a set of eigenvalues and transition amplitudes). But if N is expected (and hoped) to be very 
large (a desire a physicist has when designing a fundamental theory), the divergences can only be removed by making 
the coupling constant g vanishingly small, and fine tuned to the value of A = . On the other hand, a good 
model is expected to capture physics in a natural way, without worrying about very large numbers that by no means 
are comparable to the scale of observables of immediate interest. So, the trouble with the diverging model can be 
summarized as follows. 

Initially (in the model, by assumption; in realistic theories, motivated by experience gathered through observations 
and experiments), one is led to believe that certain Hi is a good candidate to consider. For example, in the model 
one may imagine that one measures transition rates between two states, say \'mi) and \m2), and that these rates can 
be described in first-order perturbation theory by an interaction Hamiltonian with matrix elements 

{m2\Hi\m2), {m2\Hi\mi) 
_ (mijiJ/jma), (mi|i?/|rni) 

(16) 

One is then compelled to postulate that the whole matrix of Hi has the form given in Eq. ([4]). This way of thinking 
leads to divergences, as described in the previous section, and the question what to do about them. It is clear now that 
one does not want to abandon the proposed interaction entirely since it does work for small g in the cases of interest 
(in the example, for evolution of states built from |mi) and |m2)) and exhibits appealing symmetry. In fact, this way 
of thinking is used in building theories, by extrapolation from known examples. An analogy in QFT is provided by 
the gauge symmetry [12] . Thus, one needs to conceive a general way out of the problem with divergences that are 
produced by naive extrapolation of knowledge from a small set of matrix elements to a large set of them. The large 
set is desired when one seeks a theory of presumably large range of applicability and a lot of predictive power. 

The basic idea described here (using the model) has been originally formulated in Refs. [H, The idea is to learn 
what happens when one starts with some large N and tries to reduce the value of to a smaller value, say A^i, and 
properly includes all effects due to states between N and Ni. (This is sometimes called "integrating out high-energy 
degrees of freedom.") The resulting Hamiltonian limited in energy by Ai = b^^ , will contain information about 

what has to be done in order to compensate for the presence of the arbitrarily chosen cutoff Ai . The step of reducing 
the size of a cutoff can be repeated. One can reduce Ai to A2, A2 to A3, and so on. Such steps are called RG 
transformations. When a transformation is applied K times, a chain of Hamiltonians is obtained, including H, H'^^\ 
H^^\ Eventually, two things happen. 

First, the relationship between iJ*^^) and may become universal in the sense that it no longer depends 

on all details of the initial H. For example, if contains an interaction term of the form —g'^^^y/EmEn and 

fj{K+i) contains an interaction term of the form —g^^^'^^^EmEm the relationship between 5'^' and (^(-^+1) may be 
independent of the value of the initial coupling constant g in H . Instead, a universal recursion is found for the coupling 
constant when the cutoff is changed. (Such recursion may include rescaling variables in order to compare successive 
Hamiltonians in terms of functions of dimensionless variables in a fixed range, leading to anomalous dimensions.) 
We will see how the recursion emerges in the model shortly. The model illustrates this way how the existence of a 
/3-function in QED [l^ could be understood, and explained in Ref. [l^ in the context of strong interactions. This is 
also how universality in critical phenomena could be explained in classical statistical mechanics [llj | . 

Second, after many RG steps, one obtains Hamiltonians with running cutoffs Kk- Suppose that one reduces the 
cutoff in every RG step by a factor of h. Then Kk = h^~^ may be finite even when A is sent to infinity. What is 
required is that the number of the RG steps, K , increases when one increases N for fixed N — K. Since the finite 
cutoff Kk can be chosen arbitrarily and, by construction, eigenvalues smaller than Kk for all K considered do not 
depend on at all, one eventually obtains a family of effective Hamiltonians, Hx labeled by a finite cutoff parameter 
A = h-K- Predictions that follow from Hx do not depend on A. 



Er, 



5 



A. Gaussian step 



It is time now to attempt the RG procedm'e described in the previous section in the case of our model. We will see 
how lowering A leads to increase of a coupling constant and why this increase causes trouble. 

There are essentially two interrelated tasks to accomplish. One is to establish the Hamiltonian with the large cutoff 
N, so that the resulting theory does not produce divergences (dependence of observables on TV). The other one is to 
evaluate H\ with a finite A from the well-defined initial H with a priori arbitrarily large N. The conceptual difficulty 
of the RG procedure is that the first task is accomplished in the process of trying to complete the second one, and 
in the second task, the result for H\ is used to decide how the initial H should be defined in order to make sure 
that matrix elements of H\ do not depend on A. The need for executing this process in a sequence of successive 
approximations that can converge on a structure that one is looking for, is the hardest aspect of the procedure to 
understand. 

In the simple model, one can start with the eigenvalue problem 

H\ij) = E\ij), (17) 

N 
k=M 

in which the coefRcients ipk (a wave function) satisfy the set of linear equations ([TO)) . One RG step is done by 
eliminating ipN from the remaining N — M equations. Note that the initial number of equations is iV — M + 1, 
because there is one equation with i?o = 1, in addition to N equations corresponding to positive and —M equations 
corresponding to negative powers of 6. 

The eigenvalue problem is split into one equation for the highest energy component and a set of equations for the 
remaining components: 

N 

En^N + ^ HiMnlpn = Ei/jn , (19) 
n=M 
N 

Em^m + Hiran^r. = E^^m , M <m<N -I. (20) 

n=M 



Eq. p9)) produces 

N-l 

iPn ^ [E - En - Hinn)~^ J2 ^iNn-ipk , (21) 

which is used in the remaining N — M equations to eliminate from them. In its essence, the RG step is Gaussian 
elimination of one equation in a set of linear equations. The set involves the unknown eigenvalue E. The result of 
the first step is a set of equations with M < m < N — 1, 

N-l N-1 ^ ^ 
Emipm, + Himn V'ri + T; ^7 77^ V'n = Elpm ■ (22) 

E ~ En - tiiNN 

n—M n—M 

Therefore, the new interaction "Hamiltonian" in this one-step smaller set of equations has the following matrix 
elements (M < m,n < N ~ 1): 

„(1) _ „ HiraNHiNn ,„„^ 

Himn - H,„,,, + ^-^—^ . (23) 

The word "Hamiltonian" is used in quotation marks because the matrix elements of Hj^^ depend on the unknown 

eigenvalue E. Nevertheless, Eq. (p3)) is an exact result in the model. It guarantees that H^^^ = Hq + Hf^ with the 
cutoff Ai = h^~^ has the same eigenvalue E as the Hamiltonian H with the cutoff K = h^ has. On the other hand, 
the operator ijj^^ is not fully defined before one specifies how to find the eigenvalue it depends on. Some consistency 
conditions would have to be imposed in a way that still allows us finding the initial Hamiltonian and reliably calculate 

The situation simplifies considerably if one can limit the RG procedure to eigenvalues E that are much smaller than 
a suitable finite cutoff A one wishes to reach. "Suitable finite cutoffs A" in realistic theories are cutoffs that are small 
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enough so that one can solve the eigenvalue problem for H\ on a computer. This condition puts severe constraints on 
the size of A in realistic theories, illustrated by, e.g., computational limitations of the lattice gauge theory. One has 
to execute N — n relatively complex RG steps in order to reduce the cutoff from a formally infinite A = &^ to some 
finite A and obtain ifj^ in H\ with A = 6" that is sufficiently small for a reliable computation of the eigenvalues. 
H\ still includes a large number of matrix elements, on the order of {n — AI + 1)^ times a potentially large number 
of other degrees of freedom besides the size of momentum, such as the numbers of virtual particles and their angular 
momenta, spins, colors, or flavors. 

In order to see how the simplification mentioned above emerges for small eigenvalues E, and how the simplification 
eventually ceases to be valid in the case of asymptotically free Hamiltonians, one may step back to Eq. ([221) and check 
what happens in our model. Conclusions will not be limited to the model case. 



B. Asymptotic freedom in the model Hamiltonian 



In the model, Eq. 



reads 



H 



(1) _ 



{-g\/ErnEN) {-g^/ENEn . 



g^E_ 



N 



E-E 



N 



E — En + gEff 



V E„iEn 



gE_ 



N 



(24) 
(25) 



It is clear that the interaction "Hamiltonian" Hj^'^ has the same structure of matrix elements as Hj but contains a 
new "coupling constant," say g^^\ that depends on the eigenvalue E. g^^^^ is given by the expression in the bracket in 
Eq. (Uni). 

The simplification for cutoffs much larger than E becomes obvious when one re-writes Eq. (1251) as 



^(1) - 



g^EN 



E-E 



N 



gE 



N 



= g 



1 - E/E_ 



N 



1 - .g - E/E_ 



N 



(26) 



One can neglect the ratio E/En provided that the eigenvalue E is small in comparison to the energy En = . With 
this simplification, Eq. ([26|l reads 



and implies the following recursion in further RG steps 



g 



1-5 



7(^-1) 



l-g{K-l) ' 

for as long as E^-k ^ l-^l and 1 — 5 » \E\I Efq_ji- The recursion of Eq. 

^ i-gK 



is solved by 



(27) 



(28) 



(29) 



uN-K 



and establish that the coupling 



Suppose one can solve the eigenvalue problem for H'^^^ — with Aq 

constant ^(Ao) = g'^^^ should have some value 50 in order to reproduce some measured eigenvalue _Eo ^ Ao. Eq. 
says that 



50 = 



5A 



1-^lnA/Ao 



(30) 



Eq. ([50)1 allows one to calculate that needs to stand in the initial H with the cutoff A to produce the same 
eigenvalue. The result is 



5A 



50 



l + j|LinA/Ao' 



(31) 



This means that the model is asymptotically free: the larger the cutoff A in the initial _ff, the smaller the coupling 
constant in it. The most famous analogy in QFT is found in Refs. pTi. [isj. 
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C. 



Calculation of counterterms 



Since Eq. ((3T|) has been obtained from an exact RG procedure under the two conditions that \Eo\ ^ Aq and 
1 ~ 5o ^ l-E'ol/'^O: one may expect that the replacement of g by given by Eq. pTjl provides a model H with all 
eigenvalues E that satisfy the same two conditions being also independent of A. Thus, these eigenvalues are also 
expected to be common for all Hamiltonians H\ with A > Aq, i.e., they do not depend on the finite running cutoff A 
that limits the momentum range in H\. 

The model calculation presented so far provides an example of how one takes the first step in identifying the initial 
H. The point is that the matrix elements of the initial H can be now written as 



where the letters CT are chosen as an abbreviation for the word "counterterm." The counterterm guarantees that 
small eigenvalues do not depend on the ultraviolet cutoff A. The size of the coupling constant g that was considered 
small in Section IIIIl can now be understood as actually concerning g^ in the initial, regulated Hamiltonian that 
includes counterterms; a bound state exists if g\ is sufficiently large and this means different sizes of g^ for different 
values of A. 

Note that the finite part of the counterterm, which amounts to the choice of at some value of A, is fixed by the 
condition of matching one eigenvalue with experiment. This one condition by itself (fixing one parameter to match 
one energy) would not guarantee that other eigenvalues are cutoff independent. But it does provide such guarantee 
when all matrix elements of the CT (many numbers) are calculated in the RG procedure and fixing just one value 
of the coupling constant properly determines the entire Hamiltonian matrix so that the low energy eigenvalues do 
not depend on the cutoff A (and, by construction, also do not depend on the cutoff A). One may expect that other 
matrix structures than just \JEmEn emerge in H\ with small A, but these structures must cease to depend on A when 
A — > oo once is set according to Eq. ((3T|) . 

In more detail, the argument that other eigenvalues of Elx^ that are much smaller than Aq will also be independent 
of the cutoff A once one inserts in the initial H the counterterm defined in Eq. is following. Eq. ([^5]) is valid 

and the same for eigenvalues that satisfy the conditions \E\ ^ Aq and 1 — 50 ^ |£'|/'^o- In fact, the exact way 
in which all matrix elements of H^^^^^^^^E) in the simple model depend on the initial cutoff A is determined by 
the dependence of g'^' on the sequence of ratios E/A, bE/A, 5^~"i?/A = E/X and the initial value of the 
coupling constant, i.e., after inclusion of the counterterm, just ^a- For as long as these ratios are negligible, all matrix 
elements of "Hamiltonians" H^-^^E) do not actually depend on the value of E and do not depend on the cutoff 
A. Therefore, their eigenvalues are also independent of A. However, one cannot claim that all their eigenvalues are 
the same because the number of eigenvalues depends on the size of the cutoff (the smaller the cutoff the smaller the 
number of eigenvalues). Only the eigenvalues much smaller than the smallest cutoff A in the recursion are the same 
and independent of A when calculated from all Hamiltonians Hx in which the coupling constant changes from on 
value of A to another according to the asymptotic freedom formula. 

Additional steps in evaluation of counterterms would be required in order to obtain corrections to the leading 
counterterm in Eq. (j33D due to the ratios E/Ak differing from in the RG recursion. Readers interested in such 
steps may consult Ref. 19]. It is not necessary to discuss these additional steps here in order to see the problem with 
increasing of g\ when A decreases. 

However, before proceeding to the issue of large coupling constant at low energies (in the next section), it should 
be mentioned for readers used to thinking about RG in terms of differential equations, that the reduction of one 
running cutoff on momenta to another one that is smaller, can be made in infinitesimally small steps (in the case of 
continuous momentum variables) . The corresponding changes in the Hamiltonian matrix elements are then described 
by differential equations. 



mn 



mn 




(32) 
(33) 



D. Difficulty with gx that grows at small A 



Consider Eq. (|26p which implies the recursion 



9x/biE) = gx{E) 



I- E/X 



(34) 



\^gx{E)~E/X ■ 



As discussed in the previous section, this recursion simplifies enormously when on can disregard the ratio E/X. And 
when one does, one obtains the conclusion that g\ increases when A decreases (this is equivalent to asymptotic freedom 
which says that g\ decreases when A increases). 
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The increase of g\ toward small A causes a major difficulty because the ratio E/\ is compared with 1 — g\ in the 
denominator in Eq. ([M]) . When the initial value of gx at some large A is small and gx increases when A decreases, at 
some point the effective coupling constant approaches 1. At this point, the ratio E/X cannot be neglected no matter 
how small it is. In other words, the entire procedure can no longer be based on simplifications that rely on smallness 
of the eigenvalue in comparison to cutoff. 

Not only the simplifications described in the previous section no longer apply, but also the entire procedure becomes 
unstable because the small ratio occurs in the denominator. The small denominator is reminiscent of small energy 
denominators that occur in perturbation theory and cause infrared singularities. However, the simple model shows 
that the problem is related to the strength of interactions that are capable of canceling kinetic energy terms and 
producing negative elements on the diagonal of the Hamiltonian matrix. Such situation can occur when bound states 
are formed: attractive (negative) interactions are larger than the kinetic energies. In order to control what happens 
when gx approaches 1, a different RG procedure seems required. 



In the standard RG procedure, one evaluates effective Hamiltonians Hx^ with small, running cutoffs A, as is 
illustrated in the simple model in the previous sections. One finds counterterms in the initial Hamiltonian with 
cutoff A by demanding that all matrix elements of the effective Hx are independent of A and by fixing their finite 
parts by comparison with experiment (the comparison may include conditions of symmetry [20|). Then the effective 
Hamiltonians can in principle be calculated for sufficiently small A to carry out computations of their eigenvalues 
(or other observables) using computers. Unfortunately, for feasibly small cutoffs, the effective coupling constant in 
asymptotically free theories may be so large that naive expectations based on asymptotic freedom formulae may be 
inadequate. An alternative approach is offered by the similarity renormalization group (SRG) procedure [2]| . 

In the SRG procedure, one proceeds according to similar principles as in the standard approach described in the 
previous sections. One also completes two interrelated tasks. In the one task, one finds counterterms. In the other 
task, one evaluates effective Hamiltonians. The new idea is that one does not "integrate out" any degrees of freedom. 
Instead, one changes the basis states by rotating them in the Hilbert space. The rotation is designed in such a way 
that it guarantees the new Hamiltonian Hx to have vanishing matrix elements between basis states if they differ in 
energy by more than A. This means that the matrix elements of Hx that result from the second task in the SRG 
procedure, are different from only within an energy band of width lambda along the diagonal. The algebra of the 
procedure is designed in such a way that one never encounters small energy denominators (the differential version of 
SRG procedure has the same property). Moreover, the SRG procedure allows for direct evaluation of matrix elements 
of Hx without knowing anything about eigenvalues. The reader will easily find all required details in the original 
literature. 

The SRG task of evaluating Hx with small A aims at evaluating a near-diagonal Hamiltonian matrix in such a 
way that the calculation may be carried out in perturbation theory with a small error that can decrease when the 
order of the perturbative calculation increases (perturbative approach is required in realistic theories due to their 
complexity that initially cannot be handled in any other way). One may not apply perturbation theory to complete 
diagonalization, because this would involve solving also for non-perturbative features such as bound-states. However, 
the SRG procedure can produce a Hamiltonian matrix of a small width A by "rotating out" only those interactions 
that involve energy changes larger than A and can be treated in perturbation theory. The resulting Hx must be 
diagonalized on a computer. The smaller A, the smaller the space of states required to find the spectrum in the range 
of interest. But the smaller A, the higher order of perturbative SRG evolution is required for accuracy. A compromise 
must be found, and this is not easy in QCD (see Appendix, Section fVIIIB I) . But there is also a hope that already 
first 4 orders of the calculation will be sufficient to identify key operator structures in Hx ■ 

Of course, the SRG procedure renders the same CT in the model as the one derived in Section HV C[ Eq. 
Therefore, the initial Hamiltonian with the huge cutoff A is given in the SRG also by Eq. p2|) . 



The remaining discussion is focused on how one can evaluate Hx using SRG equations in the model ll. Next 
sections will describe the outcome of these calculations and what this outcome implies regarding the role of bound 
states in evaluation of effective theories with asymptotic freedom, or limit cycles. For completeness and reader's 
convenience, a perturbative SRG scheme that can be applied in evaluation of all Poincare generators in QFT and 
used for the purpose of deriving dynamics of effective constituent quarks and gluons in QCD, is very briefly summarized 
in Appendix, Section IVlII Bl 



V. 



SIMILARITY RENORMALIZATION GROUP PROCEDURE 




(35) 
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A. Generators of the similarity transformations 

The SRG procedure is used here in its differential version. Below, prime denotes differentiation with respect to the 
parameter s = 1/A^, chosen for convenience, s increases from to cx) when A decreases from oo to 0. 
The effective Hamiltonians are generated by the formula 

H'^ = [Tx,Hx], (36) 

and Tx is called the generator. The initial condition is set at A = oo (corresponding to s = 0), 

using Eq. (|35p . Since all matrices in the model calculations, except for Hq, are functions of A, the subscript A from 
now on will be omitted everywhere, unless it is needed explicitly. The key point of the calculation is to choose the 
generator T in a way that can shed some light on the increase of an effective coupling constant at small A. 

One choice of T originates in the beautiful flow equation proposed by Wegner for Hamiltonian matrices in condensed 
matter physics [2^ . 
is reviewed in Ref. 



23, [Hi. a whole range of applications of Wegner 's and similar equations in many-particle physics 
25j . T in Wegner's equation has the form 



To = [D,H], (38) 

where D denotes the diagonal part of the matrix iJ, or Dmn = DmSmn and Vmn = Hmn (1 — (^mn)- An alternative 
equation is obtained with the choice 

Ti = [Ho,H], (39) 

which is employed in the SRG studies in nuclear physics [1^, |23, [11] . The calculation concerning gx in asymptotically 
free theories (and limit cycle) that is reviewed here [H is done with 

Tf = [G,H], (40) 

where 

G - fHo + {l-f)D. (41) 

For / = 0, one has G — D, in which the diagonal part of interactions is fully included, and T = Tq of condensed 
matter physics. For f — 1, one has G = Hq, in which no interaction effects are included, and T — Ti oi the nuclear 
studies. For intermediate values of / € [0, 1], G includes interactions to an intermediate degree, correspondingly, and 
one can inspect what happens in various cases by varying /. 

How does Eq. (|36|) work? Since T is a commutator of Hermitian matrices, it is anti-Hermitian and generates a 
unitary rotation of H , which means that traces of all powers of H are constant. In particular, 



{TrH^)' ^ Dl] +[ V |y„„M =0. (42) 




This means that the off-diagonal matrix elements decrease if diagonal matrix elements increase, and vice versa. 
Eq. p6|l implies 




, j = 4^ An(G„, - G„) I Vran P (43) 
/ mn 

= 2^(An - D^){G^ - Gn) I K«n ? (44) 

mn 

No negative terms appear on the right-hand-side when all differences D,„ — £)„ and Gm — Gn always satisfy the 
condition 



(G,„-G„)(i?„, -AO >0 



(45) 
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Using notation AiJomn = Em — En and AHjmn = Himm — -ff/rm, this condition can be rewritten, for every pair of 
diagonal elements number m and n, as 

[fAHo + (1 - /)(AHo + AHi)]{AHo + AHj) > . (46) 

Dividing by AHq > for m > n, one obtains that if u = AHj / AHq satisfies the condition 

[f+il-f){l + v)]il + v)>0, (47) 

then the sum of squares of tlie diagonal matrix elements of H increases and the off-diagonal matrix elements of H 
decrease. (By the way, in the continuum limit for Hamiltonian matrices, such as 6 ^ 1, Eq. (j47|) provides a condition 
on a derivative of the diagonal matrix elements of Hj with respect to diagonal matrix elements of Hq.) 

For / = 0, the condition (|T7)l says that {1 + v)'^ > 0. This is always true and Wegner's generator always diagonalizes 
Hamiltonian matrices because the SRG evolution stops first when all |VmnP are zero, except for elements Vmn for 
whose subscripts Dm = Dn', these may in principle stay constant unless they change due to coupling with other 
non- vanishing off-diagonal matrix elements. 

For / = 1, the condition ([T7|) reads v > —1. This means that the diagonal part of the interaction must not decrease 
faster along the diagonal than the free energy increases. Convergence of to a diagonal matrix may fail if 

for some momenta. Since Em increases monotonically with m, the lack of convergence may occur when Hjmm rapidly 
decreases with m. This happens when a negative matrix element appears on the diagonal among positive ones, leading 
to a negative eigenvalue that corresponds to a bound-state. The negative diagonal matrix elements on the diagonal 
do not guarantee that the SRG transformation stops driving off-diagonal matrix elements to zero, but it indicates 
that bound states may interfere with convergence of the SRG evolution of matrices H\ . 

For intermediate values of /, two sufficient, mutually exclusive but not necessary conditions for SRG evolution to 
bring H to the diagonal (outside regions of degeneracy mentioned earlier), are 

AHi 1 AHi , , 

In the model, it happens that these conditions can be violated when a bound state exists. The SRG evolution continues 
to bring H to the diagonal, but the effective coupling constant may become very large. This will be explained in the 
next section. 

VI. SOLUTIONS FOR Hx 

Formula (|36p produces a set of coupled nonlinear differential equations for all matrix elements of if, which means 
{N ~M + \f functions of s, 

N 

= (Gn - GkWnkVkn , (50) 

k=M 

^mn ~ ^{Gm ~ Gn){Dm — Dn)Vmn 

N 



{G„i + Gn — '^Gk)VmkVkn ■ (51) 



k=M 

This is a formidable set of equations to solve. For example, if 6 = 2, one needs 37 states to span the energy range 
between 6*^ ^ 1 keV and ^ 70 TeV, and this implies 703 functions of s for real symmetric matrices (a complex 
Hermitian Hamiltonian would imply 1369 functions). The only known way to learn precisely how the solutions to 
these equations look like with the initial condition at s = set by Eq. p7|) . is to solve them numerically 

Fortunately, the numerical analysis (it is too extensive to review here; readers interested in the numerical analysis 
need to consult the original literature) produces results that can be summarized by a very simple, qualitatively 
accurate analytic formula for H\ 
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which becomes exact in the hmit {Em + En)/X 0. This means that the matrix elements evolve with A in a coherent 
fashion and a relatively small number of simple functions of E,,n and E„ is sufficient to reasonably well describe the 
evolution of all of them. In particular, for and £'„ <C A, the exponential factor is equivalent to 1 and SRG 
evolution of the entire low-energy corner of the Hamiltonian matrix is described by just one function denoted by g\ in 
Eq. (|52|) . In analogy with the Thomson limit in QED, this function can be extracted from the lowest energy diagonal 
matrix element of the interaction term in Hx, H\mm = Em — gxEj^i. The result is that the following quantity is 
called the effective coupling constant in a Hamiltonian with a finite SRG cutoff A, 

5a = i - Hxmm/Em ■ (53) 

An alternative definition of g\, with the same result, could be based on the interaction matrix elements between two 
basis states corresponding to the eigenvalues Em and Em+i, in analogy with the example illustrated by Eq. (I16p . 
Namely, 

g\ = —H\MM+i/ \/ E]^iEm+i ■ (54) 



Asymptotically free effective interactions for different generators 



The effective coupling constant defined in Eqs. ((53)) or ([54]) . is a function of A. The observation made in Ref. f]\ is 
that the function one obtains depends on the choice of the parameter / in the generator of the SRG transformations. 
A generic example is shown in Fig. [1] In this example, 6 = 4, iV = 16, A/ — —25, and g^^ — gA ^ 4/100, where 
A = 4I6 ~ 4 • 10^. The bound-state energy, ~ 8 • 10"*^. 

The mathematical mechanism of SRG transformations by which the effective coupling constant g\ increases to huge 
values for f — 1, or stays near 1 for / = 0, or increases to intermediate values for intermediate values of /, is explained 
in detail in Ref. f]\. The physical essence of the explanation is that SRG transformations with different values of / 
place eigenvalues of H\ on the diagonal when A — > in different places. A natural location on the diagonal for an 
eigenvalue E would be a diagonal matrix element in the basis state |m) to which Hq assigns the eigenvalue Em ~ E. 

When / = 0, interactions are fully accounted for in the generator through G, see Eq. ([4T|) . The reduction of A 
brings H\ nearer its diagonal including the interaction energy that is responsible for binding. So, some diagonal 
matrix element can naturally become negative for some basis state |m) with Em \E\ = Eb- 

When f — 1, interactions are completely ignored in G, i.e., G = Hq. The diagonal matrix elements can still 
approach from above. But at some point, the monotonic increase of diagonal matrix elements along the diagonal 
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FIG. 1: Rise of the asymptotically free coupling constant g\ at small A, drawn as a function of In A/ In 6 for 6 values of /; / = 0, 
0.5, 0.75, 0.9, and 1 (the larger /, the higher curve). The rapid increase of g\ for / = 1 below In A/ In 6 ^ —8, occurs because 
A decreases below the scale of binding energy Eb- When A — > 0, (?a for / = 1 does not actually reach 00 but \EB\/b^' ~ 10^". 
For / = 0, the huge increase of gx is absent and instead gx never exceeds order 1. See the text for explanation. 
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toward large energies (the greater subscript m of a rotated basis state the greater the diagonal matrix element) is 
broken and a decrease along the diagonal (instead of increase) occurs at some state |m). At this place, the SRG 
transformation with f — 1 stops reducing the size of the nearby off diagonal matrix elements. Instead, their size 
increases and maintains the strength required for reproduction of the bound-state eigenvalue on the diagonal for some 
state \m) with Em < Eb- The shift of the eigenvalue —Eb of a fixed magnitude to a lower-energy basis state requires a 
corresponding increase of the interaction strength. For / = 1, the shifting of = —Eb toward low energies continues 
to the very end of the SRG evolution at A = and places E = —Eb at the state with the lowest allowed E^ in the 
model: the bound-state eigenvalue appears as the diagonal matrix element {M\H\\M) when A <C 6*^. As a result of 
this cumulative shift away from the natural momentum scale for the bound-state wave function, the coupling constant 
gx continues to increase in order to eventually produce —Eb on the diagonal through (1 — gx)h^^ . Thus, g\ increases 
toward EB/b^^ , which is a huge number for a large negative M. 

When / has an intermediate value between and 1, the cumulative shift stops at certain state |m), for which 
fEm -l- (1 — /)(1 — g\)Em becomes negative and reproduces the eigenvalue —Eb- This requires gx with A ~ Em to 
increase only to ^ Es/Em- 

The above scenario of how the increase of occurs, depending on the value of /, is reflected in Fig.[T] Maximum of 
the curve with / = 0.75 occurs at —9, instead of —8.5 in the case of / = 0, and the maximal value of gx for / — 0.75 
is 4, or 6 times larger than in the case with f — 0. The maximum of the curve with / = 0.9 occurs at —10, instead 
of —8.5 in the case of / = 0, and the maximal gx for / — 0.9 is 16, or times larger than in the case with / = 0, 
and so on. The case with / = 1 leads to apparently indefinite and accelerating increase that smoothly continues 
the inverse logarithmic growth that characterizes asymptotic freedom above the scale of binding. The transition 
from an asymptotically free KG behavior to the behavior dominated by bound states occurs when gx crosses 1, cf. 
Section [IVB 

In summary, the model demonstrates a possibility that an unlimited increase in the asymptotically free interaction 
at small energies is caused by omission of interactions responsible for the existence of a bound state in the generator 
of the SRG transformations. When the generator fully accounts for the interactions responsible for binding, the 
magnitude of the coupling constant never significantly exceeds 1. 

The model example is potentially important because it suggests that the increase of the coupling constant in 
asymptotically free theories that is observed in perturbation theory (corresponding to G = Hq and / = 1) may be 
not related to the phenomenon of confinement but to the phenomenon of formation of bound states. The latter 
is in principle a simpler one to handle in theory than the former. If it were indeed the case, application of the 
SRG procedure to QCD may help theorists in penetrating the range of momentum scales near A-qcd and explaining 
hadronic states without immediate necessity to solve the problem of confinement that is relevant at distances much 
larger than the size of individual hadrons (see also Appendix, Section IVIIIBp . 

B. Comment on asymptotic freedom and limit cycle 

An apparently very small alteration of the model introduced in Section [IT] leads to new ways of thinking about 
asymptotic freedom in terms of limit cycle. RG limit cycles were discovered in the context of strong interactions in 
Ref. [Tg'I, using the nowadays standard RG procedure. The SRG approach to limit cycle is based on Ref. The 
initial Hamiltonian matrix of the altered model is 

H^rnn = \/ EmEn [ (5mn - g^- ih sgn(m - Tl)] . (55) 

where i — \/— T and the new coupling constant h can be an arbitrarily small real number. If h is an arbitrary irrational 
number, the model typically exhibits chaotic RG behavior. When h = tan ^ with p an integer greater than 2, a limit 
cycle occurs. Namely, repeating the same procedure as described in the previous sections for /i = 0, one finds that 
the coupling constant h does not evolve with A at all, while the coupling constant gx that replaces in a similar 
way as gx replaced for h = Q, oscillates with A with a multiplicative period b^ . This means that has the same 
value as gx^ if Ai — {b'P)^\2 with an arbitrary integer k. The periodicity (cycle) of the coupling constant is associated 
with existence of infinitely many bound states (in the limit M — oo) whose binding energies form a geometric series 
converging on zero with quotient 1/6^. 

RG evolution of an asymptotically free interaction forms a part of a limit cycle for all values of A for which h is very 
small in comparison to ^a- In this range, h does not matter and the Hamiltonians of the model defined by Eq. (j55p 
evolve in the same way as in the case of Eq. ([55)1 (see Fig. 3 in Ref. @]). Consider the bottom- up evolution in which 
A grows and the coupling constant gx decreases as an inverse of a logarithm of A until it becomes comparable with 
h. Then, instead of ^a continuing its logarithmic falloff indefinitely, h takes over and forces ^a to change sign and 
subsequently grow in size. This continues until a new high-energy bound state is reconstructed. When this happens, 
gx switches sign to positive again and grows rapidly above 1, its magnitude depending on the size of /. The switch 
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occurs at the scale A on the order of the energy (momentum) where the new bound-state energy is located on the 
diagonal. Further increase of A produces a falloff like in asymptotic freedom again, until gx decreases again down to 
the size of h. 

The range of scales for which the cycle looks like asymptotic freedom is given by the period of the cycle, characterized 
by the factor e'^^^, which can be very large when h is very small. Such behavior of the simple model is of general 
interest because it suggests that the hierarchy problem may stem from continuing perturbative analysis for coupling 
constants g\ 1 while overlooking formation of new generation of bound states due to very small, and so far unknown, 
interactions of the type represented by the coupling constant h. But in order to see their presence, one has to use 
SRG procedure with the generator that includes interactions in G. 



VII. CONCLUSION 



The simple model study shows that the SRG procedure may be a suitable tool to handle the increase of the coupling 
constant g\ in QCD when A — > Aqqjj. If the generator of SRG transformations does not include interactions in G, 
/ = 1 in Eq. (|¥T|) . the effective coupling constant increases to very large values quickly as soon as the SRG scale 
parameter A becomes comparable with the momentum scale that characterizes formation of a bound state. In QCD, 
the corresponding momentum scale would be much larger than the scale associated with confinement because the 
size of a single hadron is much smaller than the distances at which confinement matters. If the generator of SRG 
transformations includes interactions in G, the SRG parameter A can be brought down to the momentum scale that 
characterizes bound states and the coupling constant docs not increase to large values. These model results suggest 
that the SRG procedure should be applied to QCD because it may offer help in understanding the binding mechanism 
for quarks and gluons using expansion of H\qcd in powers of an effective coupling constant and without need for 
prior understanding of confinement. Interestingly enough, the SRG procedure may also help us establish a connection 
between asymptotic freedom and limit cycle. In this respect, the model shows that in order to handle the case of 
limit cycle the generator of SRG transformations must include interactions in G. 
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VIII. APPENDIX 



The Appendix describes examples of QFT counterparts of the concepts concerning RG procedure and bound states 
that are introduced in the main text in a simple model. 



A. Potentials with 5-functions 



When one calculates corrections to the Coulomb potential in QED, one obtains a (5-function as the Uehling term in 
Hi. This term contributes a small part in the Lamb shift. The Coulomb potential, —a/r, in the Schrodinger equation 
is changed to |29| 

V{r) = ---J^/'Hr) (56) 

The correction appears suppressed by a ~ 1/137 in comparison to the Coulomb potential. In first-order perturbation 
theory, one obtains truly tiny corrections 10"'' eV for 25* states [29]). However, the second order correction 
involves multiplication of the (5-function by itself and produces infinity. The problem is not merely due to the use of 
perturbation theory instead of solving the Schrodinger equation exactly since the (5-function (or a similarly singular 
function) leads to so strong a potential that the wave function collapses unless the singular function is somehow 
replaced by a regular one. One can attempt to derive a regular expression in perturbation theory, say, by limiting 
the range of momenta in intermediate states from above by the electron mass me times the speed of light. On the 
other hand, a complete analysis should include the formation of a bound state, and bound states are not describable 
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in perturbation theory. Some form of an effective theory is necessary [sO*] (see also Ref. [3l|). In contrast with other 
approaches, the singular potentials that are obtained in QCD (or QED) in the SRG approach, are always effectively 
regulated by form factors of width A, the width playing the role of a renormalization group parameter (e.g., see 
Ref. 32]). 

Studies of (5-function potentials by physicists have a long history [H, [H, [H, [H, [13, [H, [H, HO, Slj and there 
exist mathematical textbooks on the subject [i^jISl- The simple model used in this lecture can also be derived by 
discretizing momentum on a logarithmic scale in the S'-wave Shrodinger equation for a particle moving on a plane in 
the presence of a potential proportional to a (5-function (see the original literature). 



B. SRG procedure in QFT 

Since the model discussed in this lecture appears deceptively simple, this section provides a tele grap hic overview of 
how the SRG procedure applies in QFT. The original SRG application to QCD was outlined in Ref. [4J], for calculating 
matrix elements of light-front (LF) effective H\qcd using perturbation theory. A perturbative calculus for H\qcd 
with A 3> Aqcd in terms of creation and annihilation operators, was developed in Ref. ^sj . and shown in Ref. [46| 
to be able to produce in simple scalar theories not only the Hamiltonian operator (the energy-momentum tensor 
component ) but also other generators of the Poincare algebra. This calculus produced LF QCD Hamiltonians 
with running coupling constant gx and recently led to a reasonable description of heavy-quarkonium spectra, still using 
crude and as yet unverified approximations concerning terms order gf^ [32| . The calculus is invariant with respect to 
7 independent Poincare symmetries (including 3 boosts), satisfies required cluster property 'Sl'l, and guarantees that 
the resulting Hamiltonians Hx have form factors in the interaction vertices of width A in energy, the width being the 
SRG parameter. All these features are required when one attempts to derive the parton model and spectroscopy of 
hadrons from one and the same Hamiltonian formulation of QCD. 

The point of departure is a canonical LF Hamiltonian in QFT, which requires regularization and counterterms, 

H ^ [Hean + HcT^g ■ (57) 

The canonical gluon field A in — gauge (LF in the Minkowski space-time is defined by the condition a:+ = 
-\- =Q and — means that ^ A'^ = 0) and quark field -07 are expanded into their Fourier components at 
= 0. The Fourier components correspond to canonical creation and annihilations operators for bare quarks and 

gluons (or other bare particles in other theories than QCD). These canonical operators, say qcam are related by a 

unitary transformation J7a to their counterparts for effective particles, say ^a, 

qx = UxqcanUl, (58) 

and Uco — 1. The Hamiltonian H — H(qcan) is assumed equal to Hx{qx), but the coefficients in expansion of Hx{qx) 
in powers of qx are different from the coefficients of expansion of H(qcan) in powers of qcan- The object of the 
SRG calculation are the coefficients. One actually works in the constant basis in the operator space and evaluates 
Hx = Hx{qcan)- The SRG equation is 

^Hx = [Tx,nx], (59) 

where 

Tx = U'^Ul , (60) 
and Ux corresponds to Ux- Assuming the initial condition of Eq. ([57)) . 

Uoo = [Hcan + Hcrleg : (61) 

one derives Hx from the formula 

nx = n^+ f dLo[T^,n^] (62) 

oo 

order by order in perturbation theory, eventually replacing qcan by ^a- 

The key is Tx- TYa is split into a bilinear term in qcam denoted by Tii, and the remaining terms, denoted by 7^2- 
7^2 is assumed to contain the form factors fx in vertices (the letter / used here has nothing to do with the letter / 
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used in G in Section |VA[) . Thus, 7^2 = fG2, and the form factor is defined for all operators in the same way as in the 
following example: 

g = y[123].g(l,2,3)gLi'zL29ca„3, (63) 

fg ^ y'[123]/(123)5(l,2,3)gLigL„2«can3, (64) 

/(123) = exp [-{Ml^ - Mlf/X''] . (65) 

Namely, when an operator g contains a product of m creation and n annihilation operators, the form factor operation 
inserts the same function / of the difference between invariant masses squared, Ai^^ ~ Cl2iLiPi)'^ and M^^ — 

(^^j^i kj^ , that are associated with the particle momenta using eigenvalues of tJi. The SRG generator in QFT is 
defined by the commutator ^] 

[T,g,] = [{i-Dg^r ■ (66) 

This choice guarantees that the resulting interactions are connected (cluster property) and perturbation theory for 
evaluation of Tix never contains small energy denominators, provided the function 1 — / vanishes faster than linearly 
in the difference of invariant masses (for the function adopted in Eq. (|65p . 1 — / vanishes quadratically). 



C. Questions concerning SRG and AdS/CFT in QCD 

A few questions are posed here regarding the possibility that the SRG procedure in QCD and a RG interpretation 
of AdS/CFT duality in the context of QCD, can be related. These questions reflect how little is understood about 
the SRG procedure and formation of bound states in QCD. 

SRG procedure is expected to render Hamiltonians Hx whose eigenstates in the case of LF QCD (as a part of 
the standard model) can be used to define hadrons and calculate hadronic observables. A hadron eigenstate is a 
linear combination of Fock components |n) with various values of n, where n denotes the number of operators q\ in a 
product that is used to create the state \n) from the LF vacuum state |0). The corresponding wave functions t/jux will 
depend on A while the eigenvalue and the eigenstate as a whole will not. Expansion of a hadron state into components 
with various numbers of effective constituents will be very broad in energy for A 3> Age d (corresponding to a parton 
model state for probes acting with large energy transfers) and may be dominated by the lowest constituent-number 
components for A ~ Aqcd (corresponding to the constituent model for probes acting with small energy transfers). 

The family of Hamiltonians Hx, all equivalent, can be imagined to be related to a family of complex effective 
classical Lagrangian densities £(A, x) in the Minkowski space-time, all of them representing the same theory. The 
parameter A can be considered a 5th dimension in a larger theory and one can contemplate actions of the form 

S ^ J dXd*x£{X,x). (67) 

In such action, the coupling constant gx can be considered a function of the 5-th coordinate A, and also a function 
of X. Eventually, no hadronic observable will depend on A. Nevertheless, a rich dynamics could be considered in the 
5-dimensional space, or in more dimensions, in which the family of equivalent Lagrangians labeled by A could emerge 
as describing configurations that dominate appropriate path integrals. 

Suppose the dependence on A is analogous to Polyakov's RG interpretation of the 5th (Liouville) dimension in 
string theory [5I, 'a], cf. In the SRG procedure, however, one always deals with effective quantum field theories 
in the Minkowski space. Can these effective theories explicitly include scale dependent gluonic string dynamics that 
could be "a shadow on the wall" of what happens in the bulk? According to [34| , one may expect explicit formation 
of quantum gluonic strings in LF Hamiltonian QCD in the Minkowski space when one considers eigenstates in which 
relative motion of constituents corresponds to considerably larger energies than Aqcd and components with n ^ 1 
for effective gluons are important. Could they correspond in some way to a string theory in the bulk? 

There is another line of thinking that leads to similar questions, based on the observation that the LF wave functions 
ipn\ can be used to calculate hadronic form factors using formulae similar to the ones considered by Brodsky and de 
Teramond in Refs. with additional factors that will result from using the SRG procedure. Can one connect the 

holography proposed by Brodsky and de Teramond with the A-dependent hadronic constituent distributions (densities 
obtained using sums of moduli squared of the LF wave functions ipn\), additional SRG factors, and off-energy-shell 
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old-fashioned scattering amplitudes for effective constituents with form factors of width A, all of which would result 
from HxQCD obtained from the SRG procedure? 

Such questions will not be easily answered. For example, there is a potential complication involved in QCD due 
to an infrared limit cycle [J?]. On the other hand, the model described here shows that the SRG procedure can deal 
also with limit cycles by including interactions in the generator. 
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